Search for: All records

Dust grains influence many aspects of star formation, including planet formation and the opacities for radiative transfer, chemistry, and the magnetic field via Ohmic, Hall, as well as ambipolar diffusion. The size distribution of the dust grains is the primary characteristic influencing all these aspects. Grain size increases by coagulation throughout the star formation process. In this work, we describe numerical simulations of protostellar collapse using methods described in earlier papers of this series. We compute the evolution of the grain size distribution from coagulation and the non-ideal magnetohydrodynamics effects self-consistently and at low numerical cost. We find that the coagulation efficiency is mostly affected by the time spent in high-density regions. Starting from sub-micron radii, grain sizes reach more than 100 µm in an inner protoplanetary disk that is only 1000 yr old. We also show that the growth of grains significantly affects the resistivities, while also having an indirect effect on the dynamics and angular momentum of the disk. 

Thermal ionization is a critical process at temperatures T  > 10 3 K, particularly during star formation. An increase in ionization leads to a decrease in nonideal magnetohydrodynamics (MHD) resistivities, which has a significant impact on protoplanetary disks and protostar formation. We developed an extension of the fast computational ionization method presented in our recent paper to include thermal ionization. The model can be used to inexpensively calculate the density of ions and electrons and the electric charge of each size of grains for an arbitrary size distribution. This tool should be particularly useful for the self-consistent calculation of nonideal MHD resistivities in multidimensional simulations, especially of protostellar collapse and protoplanetary disks. 

Dust grains play a major role in many astrophysical contexts. They affect the chemical, magnetic, dynamical, and optical properties of their environment, from galaxies down to the interstellar medium, star-forming regions, and protoplanetary disks. Their coagulation leads to shifts in their size distribution and ultimately to the formation of planets. However, although the coagulation process is reasonably uncomplicated to numerically implement by itself, it is difficult to couple it with multidimensional hydrodynamics numerical simulations because of its high computational cost. We propose here a simple method for tracking the coagulation of grains at far lower cost. Given an initial grain size distribution, the state of the distribution at time t is solely determined by the value of a single variable integrated along the trajectory, independently of the specific path taken by the grains. Although this method cannot account for processes other than coagulation, it is mathematically exact, fast, inexpensive, and can be used to evaluate the effect of grain coagulation in most astrophysical contexts. It is applicable to all coagulation kernels in which local physical conditions and grain properties can be separated. We also describe another method for calculating the average electric charge of grains and the density of ions and electrons in environments that are shielded from radiation fields, given the density and temperature of the gas, the cosmic-ray ionization rate, and the average mass of the ions. The equations we provide are fast to integrate numerically and can be used in multidimensional numerical simulations to self-consistently calculate on the fly the local resistivities that are required to model nonideal magnetohydrodynamics. 

Context. The degree of coupling between the gas and the magnetic field during the collapse of a core and the subsequent formation of a disk depends on the assumed dust size distribution. Aims. We study the impact of grain–grain coagulation on the evolution of magnetohydrodynamic (MHD) resistivities during the collapse of a prestellar core. Methods. We use a 1D model to follow the evolution of the dust size distribution, out-of-equilibrium ionisation state, and gas chemistry during the collapse of a prestellar core. To compute the grain–grain collisional rate, we consider models for both random and systematic, size-dependent, velocities. We include grain growth through grain–grain coagulation and ice accretion, but ignore grain fragmentation. Results. Starting with a Mathis-Rumpl-Nordsieck (MRN) size distribution (Mathis et al. 1977, ApJ, 217, 425), we find that coagulation in grain–grain collisions generated by hydrodynamical turbulence is not efficient at removing the smallest grains and, as a consequence, does not have a large effect on the evolution of the Hall and ambipolar diffusion MHD resistivities, which still drop significantly during the collapse like in models without coagulation. The inclusion of systematic velocities, possibly induced by the presence of ambipolar diffusion, increases the coagulation rate between small and large grains, removing small grains earlier in the collapse and therefore limiting the drop in the Hall and ambipolar diffusion resistivities. At intermediate densities ( n H ~ 10 8 cm −3 ), the Hall and ambipolar diffusion resistivities are found to be higher by 1 to 2 orders of magnitude in models with coagulation than in models where coagulation is ignored, and also higher than in a toy model without coagulation where all grains smaller than 0.1 μ m would have been removed in the parent cloud before the collapse. Conclusions. When grain drift velocities induced by ambipolar diffusion are included, dust coagulation happening during the collapse of a prestellar core starting from an initial MRN dust size distribution appears to be efficient enough to increase the MHD resistivities to the values necessary to strongly modify the magnetically regulated formation of a planet-forming disk. A consistent treatment of the competition between fragmentation and coagulation is, however, necessary before reaching firm conclusions.